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Abstract 

We study the properties of color-singlet Reggeon compound states in multicolor high- 
energy QCD in four dimensions. Their spectrum is governed by completely integrable 
(l+l)-dimensional effective QCD Hamiltonian whose diagonalization within the Bethe 
Ansatz leads to the Baxter equation for the Heisenberg spin magnet. We show that 
nonlinear WKB solution of the Baxter equation gives rise to the same integrable structures 
as appeared in the Seiberg-Witten solution for N = 2 SUSY QCD and in the finite-gap 
solutions of the soliton equations. We explain the origin of hyperelliptic Riemann surfaces 
out of QCD in the Regge limit and discuss the meaning of the Whitham dynamics on the 
moduli space of quantum numbers of the Reggeon compound states, QCD Pomerons and 
Odderons. 
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1. Introduction 



In the last decades a tremendous progress has been achieved in understanding of low- dimensional 
integrable field theories and statistical models. Apart from the fact that their study serves as 
a laboratory for testing a new ideas, it has been expected for a long time that low- dimensional 
integrable systems should appear in the description of the effective dynamics of quantized (3+1)- 
dimensional Yang-Mills theories (see e.g. |l]). These expectations were confirmed in 1994 by 
two examples coming from different problems: the exact calculation of low-energy effective 
action of N = 2 supersymmetric Yang-Mills theory and study of the Regge asymptotics of 
perturbative QCD || |J. It was found that well-known integrable structures emerge in both 
problems. 

In the first example ||, the Seiberg-Witten solution for the low-energy effective action of 
pure N = 2 SUSY Yang-Mills theory with the SU{2) gauge group and its generalization to 
other gauge groups || were interpreted within the framework of finite-gap solutions of the 
soliton equations and their Whitham deformations ||. In particular, the moduli space of vacua 
of N = 2 SUSY SU (N c ) Yang-Mills theory without matter was identified as a moduli space of 
hyperelliptic curve T, which amazingly enough coincides with the spectral curve of 1-dimensional 
periodic Toda chain with N c sites || [7| . Moreover, including the matter hypermultiplets in the 
adjoint and fundamental representations leads to a deformation of the Toda spectral curve into 
the spectral curves corresponding to the Calogero-Moser integrable models |8], |9| and Heisenberg 
spin chains [fL0|| , respectively. The effective energy and the BPS spectrum of excitations can 



be expressed in terms of periods of a certain meromorphic 1-differential on T, which has a 
defining property that its external derivative is a holomorphic differential on T. This differential 
was identified as a generating functional of the Whitham equations describing the adiabatic 
deformations of soliton solutions [11[ . There were different proposals to interpret the Whitham 



dynamics || |9], |T^, [13] but, despite of a lot of efforts, the connection between integrable systems 
and the Seiberg-Witten solution still has a status of an observation and the underlying dynamical 
mechanism remains mysterious. 

Our second example is related to the asymptotics of QCD at high energy and fixed transferred 



momentum, s ^> —t, the famous Regge limit [fL4}| . It has been observed many years ago (see 



e.g. |15[) that perturbative expressions for the scattering amplitudes in QCD in the Regge limit 
have a striking similarity with predictions of a low- dimensional field theory. Indeed, it was 
realized that in the Regge limit QCD should be replaced by an effective (2 + 1)— dimensional 
Reggeon field theory |l(|, which should inherit all (still unknown) quantum symmetries of 
QCD. Recently, different attempts have been undertaken to deduce its effective action out of 
QCD |L7|, [L8|. The resulting expressions become very complicated and the exact solution of 
the S— matrix of the Reggeon effective theory remains problematic. Nevertheless, there exists 
a meaningful approximation to the Reggeon effective theory, the so-called generalized leading 
logarithmic approximation (LLA) ||19H , in which the S— matrix of QCD exhibits remarkable 
properties of integrability. It turns out that in this approximation QCD is effectively described 
by the 1-dimensional quantum XXX Heisenberg magnet of spin s — J3|, [|J]. As a result, one can 
calculate the Regge asymptotics of the scattering amplitudes by applying a powerful Quantum 
Inverse Scattering Method p0 |. In particular, the Bethe Ansatz solution of perturbative QCD 
S— matrix in the generalized LLA has been developed in [0, 2T|. It is based on the Separation 



of Variables and leads to the Baxter equation for the XXX Heisenberg magnet p2| , p3[ f24]| . 

The exact solution of the Baxter equation was established in a few special cases || |2l| and 
different approaches to its general solution were proposed in |25|, |2(| |28[. One of them is 
based on the quasiclassical expansion of the Baxter equation [J2l| . It allows to obtain asymptotic 
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approximations to the solutions of the Baxter equation, which are in a good agreement with 
the results of numerical calculations. 



In the present paper we continue the study of the Baxter equation started in [p6|| . We show 
that the WKB solution of the Baxter equation, describing the Regge limit of QCD, gives rise to 
the same integrable structures as appeared in the Seiberg-Witten solution for N = 2 SUSY QCD 
and in the finite-gap solutions of the soliton equations. We explain the origin of hyperelliptic 
Riemann surfaces out of QCD in the Regge limit and discuss the meaning of Whitham dynamics 
on their moduli space. 

The paper is organized as follows. In Sect. 2 we demonstrate the relation between high-energy 
QCD in the Regge limit and XXX Heisenberg magnet. The hamiltonian of the 1-dimensional 
XXX Heisenberg magnet appears as a kernel of the Bethe-Salpeter equation for the partial waves 
of the scattering amplitudes in multi-color QCD. Its diagonalization is performed in Sect. 3 by 
applying the method of Separation of Variables p2| , f23[ , |24| , which leads to the Baxter equation. 
In Sect. 4 we develop a nonlinear WKB expansion of the Baxter equation. We describe in 
detail the classical mechanics induced by the leading term of the WKB solution and establish a 
close relation between Reggeon compound states in QCD and soliton solutions of the KP/Toda 



hierarchy p9], p0| . The quantization conditions for the Reggeon compound states and their 
interpretation in terms of the Whitham flows are discussed in Sect. 5. The concluding remarks 
are summarized in Sect. 6. 



2. High-energy QCD and Heisenberg spin magnet 

Let us consider elastic hadron-hadron scattering in QCD at high energy, s, and fixed transferred 
momentum, t, and let us replace for simplicity nonperturbative hadronic states by perturbative 
onium states built from two heavy quarks with mass M. In this case, the scattering amplitude 
can be calculated in perturbative QCD as a sum of Feynman diagrams describing the multi- 
gluon exchanges between heavy quarks in the t— channel pi] . In the Regge limit, s ^> — t , M 2 , 
the perturbative expressions for the Feynman diagrams involve large logarithmic corrections, 
a™ ln m s, which can be classified into the leading logs (n = m), next-to-leading logs (n = m + 1) 
etc. and which have to be resummed to all orders in a s . We recall, however, that one of the basic 
concepts of perturbative approach in QCD is that one considers quarks and gluons as elementary 
excitations and, assuming that their interaction is small, takes it into account as a perturbative 
expansion in a s . The necessity of resuming perturbative series indicates that the description of 
a hadron-hadron scattering in terms of "bare" quarks and gluons is not appropriate anymore in 
the Regge limit and one has to identify instead a new collective degrees of freedom, in terms of 
which the QCD dynamics will be simpler. 

In the Regge limit, due to remarkable property of the Reggeization in QCD [J31] J, reggeized 
gluons, or Reggeons, become a new collective excitations. Although the Reggeon is built from 
an infinite number of "bare" gluons it behaves as a point-like particle. It carries the color 
charge of a gluon and it has a 4-momentum k^ = (k + , k-, k±) with the longitudinal components 
k± belonging to the plane defined by the momenta of scattering hadrons and k± being 2- 
dimensional transverse momentum. The hadrons scatter each other by exchanging Reggeons in 
the t— channel and the interaction between Reggeons is described by the effective S— matrix, 
which should be obtained out of QCD in the Regge limit. 

The Reggeon S 1 — matrix has the following remarkable properties. Reggeons propagate in 
the t— channel between two hadrons and due to strong ordering of their longitudinal momenta 



the Reggeon rapidity y = In y- can be interpreted as a "time" in the t— channel Jl6| . Then 
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the Reggeon S— matrix describes the propagation of the Reggeons in the time y. Interacting 
with each other Reggeons change their 2-dimensional transverse momenta k± and color charge. 
Therefore, although Reggeons "live" in 4-dimensional Minkowski space-time, their t— channel 



evolution is described by the effective (2 + 1)— dimensional S— matrix [17| , jig] . The exact 
expression for the Reggeon S— matrix is unknown yet and one may try to evaluate it in different 
approximations. In what follows we will study the Reggeon S— matrix in the generalized LLA 



[19|j . In this approximation one preserves the unitarity of the S— matrix in the direct (s— , t- 
and u—) channels but not in the subchannels. 



2.1. Generalized LLA 

The hadron-hadron (or onium-onium) scattering amplitude is given in the generalized LLA by 
the sum of effective Reggeon Feynman diagrams fl9j in which, first, the number of Reggeons 
propagating in the t— channel is preserved (no creation/annihilation of Reggeons) and, second, 
only two Reggeons could interact with each other at the same moment of "time" y. Therefore, 
according to the number of Reggeons, N, the scattering amplitude can be decomposed as 

oo 
N=2 

where the A^— th term describes the evolution in the t— channel of the system of A" Reggeons with 
a pair-wise interaction. The amplitude AN(s,t) satisfies the Bartels-Kwiecinski-Praszalowicz 
equation , which can be interpreted as a Bethe-Salpeter equation for the A" Reggeon scatter- 
ing amplitude. Its solutions define the color-singlet A" Reggeon compound states, perturbative 
QCD Pomerons and Odderons, whose contribution to the scattering amplitudes takes the stan- 
dard Regge form 

A N (s,t)=isJ2$ } (t)$ } (t)s E ^, 
M 

where indices A and B refer to the scattered hadrons. Here, -Ejv,{<j} is the energy of the A" 
Reggeon compound state. It is defined as an eigenvalue of the A" Reggeon hamiltonian 

Ti-N\XN,{q}) = E N,{q}\XN,{q}) , (2-1) 

with {q} being some set of quantum numbers parameterizing all possible solutions. The residue 
functions (3^ and j3g measure the overlapping between hadronic wave functions and the wave 
functions of the A" Reggeon compound states 

(3 { A q} (t) = (A\ XN , {q} ). 

The pair- wise Reggeon hamiltonian Ti^ acts in Q2.1|) on 2-dimensional transverse momenta of 
A" Reggeons and the relation ( |2.1| ) has the form of a (2+1) -dimensional Schrodinger equation. 
Let us change the representation from 2-dimensional transverse Reggeon momenta space k± to 
2-dimensional impact parameter space b± = (£, Q and define the holomorphic and antiholomor- 
phic complex coordinates of Reggeons as 

Zj = Zj + iQ, Zj = Zj-itj, j=l,2,...,N. (2.2) 

Remarkable property of the Reggeon hamiltonian TC^ is that being expressed in terms of holo- 
morphic and antiholomorphic Reggeon coordinates, it splits in the multi-color limit, a s N c = 
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fixed and N c — > oo, into the sum of holomorphic and antiholomorphic hamiltonians |32] [] 



H N = ~ff(H N + H N ) + 0(N; 2 ) , [H N , H N ] = 0, 

where Hn and Hn act on holomorphic and antiholomorphic Reggeon coordinates, respectively. 
Since the hamiltonians Hn and Hn commute with each other, the original (2+l)-dimensional 
Schrodinger equation fl2.1| ) can be replaced by the system of holomorphic and antiholomorphic 
(l+l)-dimensional Schrodinger equations 

H N \<f N ,{q}) = £ N,{q}\ i PN,{q}) , H N \Tf Nt{q} ) = S N ,{q} \^N,{q}) > ( 2 -3) 

which define the spectrum of the Reggeon compound states as follows 

E N,{q} = ^^( £ N,{q} + £ N,{q}), XN,{q}{z,z) = <p N!{q} (z)p NAq} (z) , (2.4) 

where z and z denote the full set of N Reggeon holomorphic and antiholomorphic coordinates. 
Thus, in the generalized LLA perturbative (3+l)-dimensional multi-color QCD is effectively 
replaced by the system of (1+1) -dimensional Schrodinger equations ( |2.3|) . The properties of 
two equations in ( |2.3| ) are similar and in what follows we will study only the holomorphic 
Schrodinger equation in (|2.3j ) for fixed number of Reggeons N. 



2.2. Integrability 

The holomorphic hamiltonian Hn describes the nearest-neighbour interaction between N Regge- 
ons with holomorphic coordinates z\, zn and periodic boundary conditions Zk+N = Zk 

N+l 

H N =J2 H jd+1 (2.5) 

3=1 

and Hn,n+i = Hn,i- Here, the interaction hamiltonian between two Reggeons, Hjj+i, is given 
by the BFKL kernel [|l|, H 

H 12 = H(z u z 2 ) = -^(-Ju) - V(l + J12) + 2^(1) , (2.6) 

where ip(x) = dln JJ x ^ and the operator J 12 is defined as an operator solution of the equation 

Ji2(l + Jn) = -(z 1 -z 2 ) 2 d 1 d 2 

with dk = d/dzk- We stress that the operator (|2.6f) was originally derived |3l| from the analysis 
of Feynman diagrams contributing to the Regge asymptotics of the onium-onium scattering 
amplitude in perturbative QCD. Therefore, it was quite unexpected to find || |J that it coincides 
with the hamiltonian of 1-dimensional XXX Heisenberg magnet of spin s = corresponding to 
the principal series representation of the SX(2,(D) group. The number of sites of the spin chain 
is equal to the number of Reggeons. 

Among other things this equivalence means that the holomorphic Schrodinger equation ( p.3| ) 
possesses the family of hidden mutually commuting conserved charges |3]. [| 

[q k ,q 3 } = [q k ,H N )=0 (2.7) 



1 For N = 2 and N = 3 this relation is exact even for finite N r . 
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and their number is large enough for the system ( |2.3|) to be completely integrable. To identify 
the conservation laws for the system of iV interacting Reggeons one applies the Quantum Inverse 
Scattering Method [21J and uses well-known constructions for the Heisenberg magnet [[34], |35|| . 
Namely, to each Reggeon we assign the (auxiliary) Lax operator 

( A + iSl iSZ \ 



iSt A - iSf 



A-l + i )®(z k ,-l)d k , (2. 



k llD k ) \z k 

where A is an arbitrary complex parameter and S k are spin s = generators of the SL(2) group 
acting on the holomorphic coordinates of the k— th Reggeon 

Sf = z k d k , S k = z\d k , S k = —d k . 

Then, we construct the monodromy matrix 



T(A) - Ljv(A)Lat-i(A)...Li(A) - ^ ^j^j ^j^j 



(2.9) 



where the operators A, B, C and D act on the holomorphic coordinates of A" Reggeons and 
satisfy the Yang-Baxter equation |35|, |34j| . Finally, we obtain the transfer matrix as 



A(A) = trT(A) = A(X) + D(X) (2.10) 

and verify, using ( |2.8| ) and ( j2.9| ), that A(A) is a polynomial of the degree iV in A 

A(A) = 2X N + q 2 X N ~ 2 + ... + q N , (2.11) 

with the coefficients q k being operators acting on holomorphic coordinates of the Reggeons. It 
immediately follows from the Yang-Baxter equation for the monodromy matrix T(A) that the 
operators q 2 , qN satisfy the relations Q2.71) 0- 

Thus, we identify the operators q 2 , qN as — 1 conserved charges for the system of A" 
interacting Reggeons. To match the number of conservation laws with the number of Reggeons, 
one needs one more conserved charge. It is easy to see that the remaining A^— th charge is 
associated with the center-of-mass motion of the A^ Reggeon compound state and it is equal to 
the total Reggeon momentum 

P = 7ri + 7r 2 + ... + ir N = i(S7 + 57 + ... + S N ) 

with TTj = —idj being the holomorphic component of 2-dimensional transverse momentum of 
the j— th Reggeon. The explicit expressions for the operators q k can be obtained from (|2.8|) , 
(PD and dTJ) as 1, | 

q k = ^2 i z jij2 z j2j3--- z jkji^ji^h---^jk (2-12) 

A r >ii>i2>->j fc >i 

with Zj k = Zj — z k . They can be interpreted as higher Casimir operators of the SL(2) group and 
their appearance can be traced back to the invariance of the Reggeon hamiltonian (|2.5| ) under 
the Mobius transformations 

az k + b 

Zk^ —:, 2.13 

cz k + a 

2 Notice that in contrast with the spin— | Heisenberg magnet and 1-dimensional Toda chain |^|, the hamil- 
tonian Hn does not enter into the expansion of the auxiliary transfer matrix (2.10). To obtain Hn one has to 
construct the fundamental transfer matrix. 
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where ac — hd = 1. In particular, 



92= £ z%d 3 d k = -h(h-l) (2.14) 

N>j>k>i 

is the quadratic Casimir operator and its eigenvalue h defines the conformal weight of the 
holomorphic wave function ipn(zi, zn) of the N Reggeon compound state (|2.4|) . The Reggeon 
wave function belongs to the principal series representation of the 5L(2,(D) group and, as a result, 
its conformal weight is quantized 

1 -\- m 

h = \-iu, , m = Z, z/ = R. (2.15) 



Here, integer m defines the Lorentz spin of the Reggeon state ( |2.4|) , corresponding to the rota- 
tions in the 2-dimensional impact parameter space, xn e ima XN a s z j — * e ta Zj and Zj — > e~ ia Zj. 



The quantization conditions for the remaining charges q^, q^ are much more involved |26 
and their interpretation in terms of the properties of Reggeon states will be given in Sect. 5. 

Once we identified the complete set of conservation laws, the Reggeon hamiltonian becomes a 
rather complicated function of the conserved charges 0, Hn = H^(q2, gjv), and one can replace 
the original Schrodinger equation ( |2.3j ) by a simpler problem of simultaneous diagonalization 
of the operators P, q 2 , qN- Their eigenvalues form a complete set of quantum numbers 
parameterizing the N Reggeon compound states, (pN,{q}{ z ii •••> z n)- 



3. Separation of variables 

Diagonalization of the conserved charges can be performed using the Separation of Variables 
(SoV) |22|, |23|, P4"| . The operators P, q 2 , qN depend on the coordinates and momenta, z k 
and 7Tfc = — idk, respectively, of N Reggeons and their diagonalization is reduced to solving 
of a complicated system of N coupled partial differential equations for the holomorphic wave 
function of the N Reggeon state. Instead of dealing with this system we perform a unitary 
transformation 

( Zj , nj) - ( Xj = U^ Zj U, Pj = U^U) , U = U({z k , n k }) (3.1) 

in order to replace the original set of Reggeon coordinates and momenta, (zj,iTj), by a new 
canonical set of separated variables 

[xj,x k ] = [Pj,Pk] = , [pj, x k \ = -i8 jk , (3.2) 

in terms of which the following relations hold 

$(xj,Pj\q2,-,qN)\<PN,{ q }) = 0, j = 1,2,..., N- 1, (3.3) 

(PN - P)\<PN,{q}) = 0. 

Here, the number of pairs (xj,pj) is equal to the number of Reggeons, N. One of the pairs, 
(xn,Pn), describes the center-of-mass motion of the iV Reggeon state 

TV N 

X N = J-rY, Z j , PN = Y, n i' ( 3 ' 4 ) 

JV 3=1 3=1 



3 Due to invariance of the Reggeon hamiltonian under Mobius transformations (2.13), Hn does not depend 
on the total momentum P. 
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while the definition of the remaining separated coordinates and the functions $ will be given 
below. In relation ( |3.3D , the operators are ordered inside the function $ as they are enlisted and 
the conserved charges q% can be replaced by their eigenvalues corresponding to the N Reggeon 
compound state <fN,{g}- 

The remarkable property of the separated coordinates is that in the a;— representation, pj = 
—id/dxj, the relations ( |3.3|) define the system of N separated differential equations for the wave 
function of the N Reggeon state fN,{q] = <PN,{q}(xi, ■■■,xn) and its solution takes the following 
factorized form 

<PN,{q}{xi, —,xn) = Q(x 1 )Q(x 2 )-..Q(x N ^ 1 ) exp(iPx N ) , (3.5) 
where P is the total momentum of the state and a function Q(x) satisfies the Baxter equation 

$(x, -id x ; q 2 , q N )Q(x) = . (3.6) 

To obtain the wave function in terms of the holomorphic Reggeon coordinates one has to go 
back from x— to z— representation in ( |3.1| ) by applying the unitary transformation U. The 
resulting expression can be found in ||. 



3.1. Baxter equation 

Let us construct the separated variables for the system of N Reggeons. Before doing this we 
would like to notice that the SoV method has been extensively developed in application to 
1-dimensional classical integrable systems. In this it was shown in many examples |p9|j . 

poles of a properly normalized Baker- Akhiezer function provide the set of separated coordinates. 
Having expressions for the separated variables in the classical system, one may try to generalize 
them to the corresponding quantum integrable model. However, one of the peculiar features of 
the holomorphic Schrodinger equation for the N Reggeon compound states is that it describes 
1-dimensional essentially quantum integrable system and one has to determine the separated 
variables without knowing their classical analogs. 

Constructing the separated variables for the N Reggeon state we follow the approach de- 
veloped by Sklyanin [^fl . Namely, we define the separated coordinates x±, xn-i as operator 
zeros of the operator B(X) entering the expression ( |2.9| ) for the monodromy matrix. According 
to (|2.8|) and (|2.9|) , B(X) is a polynomial of the degree N — 1 in A and it can be represented as 

B(X) = p N (X - Xi)(A - x 2 )...(A - xjv_i) , 



where p^ was defined in ( |3.4j) and the ordering of operators is unessential due to ( |3.2|) . The 
definition of the momenta pj conjugated to the coordinates Xj takes a few steps. Let us substitute 



A = Xj into ( |2.9|) . Since the operators A(X), C(X) and D(X) do not commute with Xj we use 
the prescription of "substitution from the left", A i— > Xj. This means ||24|1 , that for an arbitrary 
operator W(X) = J2k ^ k Wk one first pulls out all powers of A to the left and then replaces A by 
the operator Xj as W(X ^ Xj) = Y,k x jWk- For A t— > Xj the monodromy matrix ( p.9|) takes a 
triangle form with the diagonal elements given by 

A(X ^ xj) = xfoo+ , D(X ^ Xj ) = xfooj , (3.7) 

where notation was introduced for the operators uj^ and toj. 

Let us forget for a moment that we are dealing with operators and replace all operators by 
the corresponding classical functions. Then, one could identify the diagonal elements x^ujf of 
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the monodromy matrix T(A) as its eigenvalues and write down the characteristic equation for 

cj ± as 

det(T(x) - x n uj) = x 2N uo 2 - x N u trT(ar) + det T(x) = . 

Using the definition of the transfer matrix, fl2.1U|) , and that classically det T(A) = JlfeLi det L k (X) 
= X 2N we obtain the relation 

uj + - = x~ n K{x)=2 + % + ... + %, (3.8) 

UJ x z X JV 

which can be interpreted as the level surface of the commuting integrals of motion q k , or spectral 
curve of the monodromy operator for some classical Liouville integrable system |?9] . As we will 
show in the next section, this system arises as a quasiclassical limit of the N Reggeon state. 

Let us take now into account noncommutativity of operators and establish the quantum 
analog of the spectral curve ( |3.8|) . Using the Yang-Baxter equation for the monodromy matrix 
one can verify the following relations for the operators uft 



[uf, uj^\ = [uf, wjf] = , uj+ujj = 1 , ujfx k = (x k ± iS jk )uf . 

These relations suggest that the operators uf act on the Reggeon wave function in the a:— repre- 
sentation as shift operators and they allow us to define the momentum pj conjugated to Xj in 
the separated variables as follows 

uja = uf = — = ± exp(-pj) , u k Q(xj) = ±Q(xj + i5 jk ) . (3.9) 

We notice that the canonical commutation relations ( |3.2|) do not allow to fix a sign ambiguity 
in the definition of the momentum. 

To establish the relations ( |3.3| ) in the separated variables (xj,pj), we substitute A i— »■ Xj into 
the definition ( [2.10| ) of the transfer matrix and use ( p.7| ) to get 

Uj + — = x- JV A(x j ) = 2 + xfq 2 + ... + xj N q N , (3.10) 

UJj 

where j = 1, N — 1 and the ordering of the operators Xj and q k is important since they do 
not commute with each other. The operator identity ( P-lOj ) can be considered as a quantum 
analog of the spectral curve (|3.8| ). It relates a pair of the separated variables (xj,pj) to the set 
of the conserved charges and being applied to the wave function of the N Reggeon compound 
state it leads to (|3.3|), provided that 

^(x j ,p j ;q 2 ,...,q N ) = x~ N A(xj) -2coshp j . (3.11) 

Here, we used ( |3.9|) with a plus sign and took into account the possibility to change a sign by 
introducing an ambiguity pj — > pj + in in the definition of the momentum^- Substituting (ft.lip 
into ( p.6| ) we find that the wave function Q(xj) satisfies the finite-difference equation 

xj N A{x j )Q{x j ) = Q{x j + i) + Q{x j -i), j = 1,2,..., N- 1, (3.12) 

which has the same form as the equation for eigenvalues of the Baxter Q— operator P? |. In 
this equation, Xj denotes the eigenvalue of the corresponding coordinate operator. Notice that 

4 Similar to the situation in the Toda chain p^ j. the sign ambiguity can be fixed by requiring the momentum 
Pj to have real eigenvalues. 
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the operators qk entering into the definition of the transfer matrix (|2.11|) have been replaced 
in ( |3.12j ) by their eigenvalues corresponding to the N Reggeon compound state. As for any 
quantum system, the possible values of qk are constrained by the quantization conditions which 
will be discussed below. 

Having solved the Baxter equation ( |3.12| ), one can obtain the wave function of the N Reggeon 
compound state ( |3.5|) . Moreover, the same function Q(x) controls the dependence of the holo- 
morphic Reggeon hamiltonian on the conserved charges and for a given set of quantum numbers 
q2, qN the holomorphic energy of the N Reggeon state, defined in ( |2.3|) , can be evaluated as 



>Qn) 



d Q(x — i) 

l ~r ln 7^7 ^ 

ax Q[x + %) 



(3.13) 



The Baxter equation ( |3.12| ) alone does not allow us to find its solution Q(x). One has 
to specify additionally the spectrum of the operators x%, xn and provide the appropriate 
boundary conditions for the function Q(x). It is clear that the possible eigenvalues of the 
coordinate operators Xk depend on the class of functions Q(x) on which they act. For example, 
the spectrum of Xj is discrete for the spin— 1 Heisenberg magnet, while for the Toda chain the 
eigenvalues of Xj take continuous real values |22|, ^3], p4}| . To find the spectrum of the operators 
Xj for the Reggeon states we have to impose additional conditions on the functions Q(x). To 
this end we notice that the solution of the Baxter equation Q(x) defines the wave function of 
the Reggeon compound state ( |3~5| ) and the above conditions should follow from the requirement 
for its norm to be finite. Although the explicit expression for the norm is known f|], |j, the 
general form of the resulting constraints on the functions Q(x) was not found yet, except of the 
subclass of the so-called polynomial solutions of the Baxter equation [p6|| , corresponding to the 
special values of quantized conformal weight fl2.15| ). 



3.2. Polynomial solutions of the Baxter equation 



The subclass of polynomial solutions covers all normalizable solutions of the Baxter equation 
corresponding to the integer positive conformal weight of the N Reggeon compound states 



h> N . 



(3.14) 



where one put v = and m > 2N — 1 to be odd in ( |2.15| ). In this case, as was shown in pfifl . 
the normalizable solutions Q(x) of the Baxter equation become polynomials of the degree h in 
x with N time degenerate root x = 



h-N 



Q(x) 



n (* - 



(3.15) 



Each solution is characterized by the set of nonvanishing roots {Xj}- Having the explicit expres- 
sions for the roots, one can easily evaluate from ( |3.12|) and (|3.13|) the quantum numbers qk and 
the corresponding energy of the iV Reggeon states. It follows from the Baxter equation ( p,12[ ) 
that Xj satisfy the Bethe equations for the XXX magnet of spin s = |21] . Their study reveals 
the following interesting properties of the polynomial solutions of the Baxter equation plj ^| . 

For any given integer h in ( |3.14j) , the space of polynomial solutions is finite-dimensional. 
The possible values of roots A^, as well as the values of quantized qk and the energy e^, turn 
out to be real and simple 



Im A = Im q^ = Im = 



(3.16) 
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and they can be parameterized by integers ni, ri/v-2 



qk = Qkihint, ...,n N - 2 ) , s N = s N (h;nx, ...,n N - 2 ) (3.17) 

such that 

N-2 

ni, ...,njv-2 > 0, ^2n k <h-N. 

k=l 

The eigenvalues of the coordinate operators X\, £jv-i are also real for the polynomial solutions 
and they occupy N — 1 compact nonoverlapping intervals on the real axis 

ImXj = 0, Xj G [c 2 j-i, 02j] • (3.18) 

The positions of the intervals, o^, depend on the values of quantized q^, qN and they satisfy 
the following relation 

A 2 (a fc )=4af. (3.19) 

The possible values of roots belong to the same intervals ( |3.18| ) and their distribution on the 
real axis was found in the limit of large integer conformal weight, h ^> 1. Moreover, in the large 
h limit the following scaling holds 

Xj = 0(h) , q 2 = 0(h 2 ) , ... , q N = 0(h N ) , e N = 0(\nh) , (3.20) 

which allows to develop the asymptotic expansion of q^ and En in inverse powers of the conformal 
weight 

oo oo 

Qk = h k J2 Qk ( n ) h ~ l > £n = -2N In h + ]T eg (n) /i"' . (3.21) 
z=o «=0 

where k = 3, and the coefficients qf' and e$ depend on the integers n^. For N = 2 

Reggeon states all coefficients are known exactly and for iV = 3 Reggeon states they were 
calculated up to 0(h~ 8 ) order |2T| 



Comparing the asymptotic expansions (|3.21| ) for N = 2 and > 3 Reggeon states one can 
discover two important differences [26]. First, the coefficients e$ grow as factorials to higher 



orders in 1/h. The asymptotic series for the energy of = 2 Reggeon states is Borel summable, 
but for N > 3 the properties of the same series are drastically changed and it becomes non Borel 
summable. Second, for A^ > 3 Reggeon states one has to study additionally the properties of 
the "higher" quantum numbers q^, q^. One finds that for A^ > 3 Baxter equation among all 
possible polynomial solutions there are a few singular solutions, corresponding to the situation, 
when either one of the intervals ( ft.lSQ is shrinking into a point, or two intervals merge, <Jj-i = <7j. 
On the quantum moduli space of q^, qN the singular solutions are described by the "critical" 
(N — 3)— dimensional hypersurface 

EtfCgg* ...,(#*) = 0. (3-22) 
For example, for A^ = 3 Reggeon states the critical values of g 3 are just three points given by 

g 3 crit = and ± J=(l + (^aT 1 )) . (3.23) 

As we will show in Sect. 4, both properties are closely related to each other and they can be 
simply understood using the properties of hyperelliptic Riemann surfaces determined by the 
complex curve ([3.81) . 



11 



The polynomial solutions of the Baxter equation are closely related to the different systems 
of orthogonal polynomials 



21 



For N = 2 the exact solution of the Baxter equation was 
identified as a continuous Hahn symmetric polynomial. For N > 3 more complicated systems of 
orthogonal polynomials arise and their properties can be studied using the approach developed 
in IBS 



3.3. Quasiclassical limit 



The derivation of the Baxter equation in Sect. 3.1 was based on the close relation between 
Heisenberg spin magnet and the system of N interacting Reggeons in multi-color QCD. It 
is not surprising therefore that similar equations appear in different 1-dimensional quantum 
integrable systems. Namely, applying the SoV to the periodic Toda chain with N sites one 
obtains the Baxter equation in the form |23|, |24], |39| 



AToda(^i) Q{xj) = Q{xj + iti) + Q(xj -ih) , j = 1, 2, N - 1 



(3.24) 



where Axoda is the corresponding transfer matrix and h is the Planck constant. Comparing 
( |3.12| ) and ( |3.24| ) we notice a trivial but important difference between two models - the Baxter 
equation for the N Reggeon state, as well as the canonical commutation relations (|3.2j) , do not 
involve the Planck constant. Indeed, the Planck constant enters into the definition of the XXX 
Heisenberg magnet as a coefficient in front of the second term in the Lax operator ( |2.8|) . The 
Reggeon Lax operator (|2.8| ) and, as a consequence, the Baxter equation (|3.12j ) correspond to 
the special case h = 1. 

Thus, the system of TV Reggeons does not have any natural small parameter, which would 
allow us to perform the quasiclassical limit, Q(x) = exp(|S'(x)) as h — > 0, and obtain its 
classical integrable analog. Nevertheless, there is a simple way of introducing a small parameter 
7] into the Baxter equation similar to the one used in asymptotic solutions of the Painleve type 
I equation f4"0fl . Let us rescale in ( p,12|) the collective coordinate as x — > x/r] and introduce 
notations 

f(x) = r ] h Q(x/r ] ), q k = q kV h . (3.25) 
Then, one could identically rewrite the Baxter equation ( |3.12| ) as 



f(x + irj) + f(x -ir}) = x N A(x, {q k })f(x) 



(3.26) 



where the transfer matrix is given by ( 2.11 ) with q k replaced by q k defined in ( ^.25| ). Following 
|(J, let us look for the solutions of the Baxter equation in the form 



f{x) = exp i-S(x, {q k }) 



S(x, {q k }) = S (x) + TjStix) + 



(3.27) 



where each term is assumed to be uniformly bounded and the expansion to be convergent. 

The asymptotic expansion ( |3.27| ) of the solutions of the Baxter equation naturally appears 
in the limit ?7 — ^ and q k = fixed. According to ( |3.25| ) and ( |2.14|) , this limit corresponds to the 
large values of the conformal weight, h ~ l/rj — > oo, and it was studied in J^J. Using (|3.15|) , 
( |3.25| ) and (|3.27|) we can express S in terms of the roots Xj as 



S (x; {qk}) = -if] m ( x ~ ^jV) 



(3.28) 
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where Xh-N+i = ••• = A^ = 0. Then, for h ~ 1/rj — > oo the property of roots ( |3.20| ) implies the 
scaling A^r] = 0{rf) and S^x; {<?£;}) = 0{rf), which leads to the asymptotic expansion Q3.27|) . 



Substituting the ansatz ( |3.27| ) into ( |3.26| ) and expanding the both sides of the Baxter equa- 
tion ( |3.26 ) to order i] one obtains [|39| 



2cos1iS , q(x) = x A(x; {<?&}) , S[(x) = -Sq(x) coth.S' (x) , (3.29) 

where prime denotes a derivative with respect to x. Then, the solution of the Baxter equation 
( |3.26|) is given by 

f(x) = — — L exp (l f dS (x) + 0( V )) . (3.30) 
[smh S (x)\ 1 \V J j 

Comparing ( |3.26| ) with the Baxter equation for the Toda chain ( |3.24j ) we observe that the 
parameter r\ plays a role of the Planck constant. The ansatz (|3.30|) takes the form of the WKB 



expansion with the leading term Sq(x) being the "Reggeon classical action". Let us study the 
classical mechanics of the N Reggeon system governed by the action Sq(x). 



4. Quasiclassical Reggeon states as KP/Toda solitons 

Let us introduce notation for the function 

uj(x) = exp(S , Q(x)) (4.1) 
and rewrite the first equation in Q3.29 ) as 



uj + - = x~ N A(x;{q k }). (4.2) 

00 



This expression is the quasiclassical analog of the operator relation ( |3.10p . We notice that 



the variable x appears in ( 3.26 ) and (|4.2j) as a collective coordinate of the Reggeons and for 
polynomial solutions it takes values inside one of the N — 1 compact intervals on the real axis 
( |3.18|) . The classical analog of the Reggeon momentum in the separated coordinates, [x,p] = it], 
can be obtained from ( |3.9| ) as f] 

p = ln|w(a:)|, (4.3) 



where u takes positive and negative values in ( |4.2j ) and the modulus is needed for momentum 
p to be real. Using (|4.2|) and taking into account that uj is real on the classical trajectories, we 
obtain the relations which define the intervals of the classical motion of Reggeons as 

A 2 (x,; {q k }) > Axf , j = 1, 2, N - 1 , (4.4) 

where index refers to the j—th allowed band. This relation is in agreement with ( |3.19| ). 

5 The same quantity can be interpreted E| as a quasimomentum for the auxiliary first-order difference linear 
problem corresponding to the Reggeon Lax operator, L n tp n (E) = Eij) n +i (E) , with the periodicity condition 
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4.1. Spectral curve 



Let us continue the relation (|4.2| ) into the complex domain. For r\ — 1 it coincides with the 
spectral curve ( |3.8| ) and replacing 



N 



y = x [LO- 



CO 



one can rewrite it in the conventional form as 



N 



y 



A 2 (x) - Ax 



2N 



Aq 2 P 2N _ 2 (x) 



(4.5) 



(4.6) 



where P 2 n- 2 (x) = x + ... is a polynomial with the coefficients depending on the quantum 
numbers qk of the N Reggeon state. The curve (|4.6|) determines the hyperelliptic Riemann 
surface Y N equipped with a meromorphic differential dS defined in (|3.29| ). The genus of the 
Riemann surface, g = N — 2, depends on the number of Reggeons inside the compound state. 
For N = 2 Reggeon state, the well-known BFKL Pomeron pi|, T 2 is a sphere, 



2 • 



y 2 = q 2 (Ax 2 + q 2 ) 



while for N = 3 Reggeon state, the QCD Odderon @, T 3 is a torus, 

r 3 : V 2 = (fax + 93) ( 4 ^ 3 + Q2X + q 3 ) ■ 



(4.7) 



(4. 



The branch points of the complex curve, y = 0, correspond to u — ±1 (or p = 0) in ( f4.5|) . They 
coincide with the turning points of the classical Reggeon trajectories and define the boundaries 



of the allowed bands, x — (jj, in (|3.19| ). Then, the Riemann surface admits a representation 



in the form of two sheets of the complex x— plane glued together along the cuts running between 
the turning points &2j-i and a 2 j (see Fig. 1). One of the sheets will be called the upper and 
the other one the lower sheet. Each point on the Riemann surface Tn can be parameterized as 
Q = (x, ±), where different signs correspond to the upper and lower sheets. 




N-l 



Figure 1: The definition of the canonical basis of cycles on the Riemann surface T^. The dotted 
line denotes the part of the (3— cycles belonging to the lower sheet. The a— cycles correspond 
to the classical motion of Reggeons over the allowed bands. 

In the center-of-mass frame, the classical trajectories of the Reggeons correspond to the 
cycles over N — 1 allowed bands on the Riemann surface^, which we denote as with k — 1, 
N — l. By the definition, for each point Q on the cycle the corresponding Reggeon coordinate 

6 For polynomial solutions of the Baxter equation the number of allowed bands should match the number of 
degrees of freedom minus the motion of center-of-mass with the coordinate xn- 
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and momentum, Xj and pj, take real values. It is easy to see from Fig. 1 that the sum of all 
a— cycles is homologous to zero 

ai + a 2 + ... + a N -i ?z (4.9) 

and one can choose the first g = N — 2 a— cycles to construct the canonical basis of cycles on 
Tat: oti, o. AT— 2) Pi, Pn-2 with the intersection matrix 

aj o a k = Pj o f3 k = , a>j o (3 k = 5 jk . 

The definition of the P— cycles is shown in Fig. 1. Obviously, the introduction of the P— cycles 
has a meaning only for N > 3 Reggeon states. 

At this point we notice the important difference between N = 2 and N > 3 Reggeon states 
- the appearance of the moduli for the complex curve for g = N — 2 > 1. To identify the 
moduli we observe that the curve (|4.6|) is invariant under transformation 

x^Xx, q k ^X k q k} y^X N y, (4.10) 

which is induced by rescaling of the parameter rj in ( [3.25 ). Let us introduce the moduli of the 
hyperelliptic curve Tn for iV > 3 as 

u k - 2 = q k /(-q2) k/2 = g fc /(-g 2 ) fc/2 , k = 3, N , (4.11) 

where we used (|3.25 ) to eliminate the dependence on arbitrary parameter rj. Thus defined 
moduli are invariant under ( 4.10|) and depend on the quantum numbers of the N Reggeon 
states. It is also convenient to introduce the following variable 

T = (-g 2 ) 1/2 = v(~q2) 1/2 ■ 

According to ( |3.16|) , for polynomial solutions of the Baxter equation, u\, an d T take 
real quantized values, such that 

<? 2 = -T 2 , q k = T h u k „ 2 , (4.12) 

with k = 3, N. 

The relation ( |4.1| ) defines the following meromorphic 1-differential on the hyperelliptic curve 

, i / x du NA(x) - xA'(x) , .„ 

dS = dx\nuj{x) = -x — = ^ y -^dx. (4.13) 

lo y 

At the vicinity of two infinities on the upper and lower sheets of the Riemann surface, P±, the 
differential behaves as 

,„ X~*P± ■rj.dX 

x 



and we identify dSo as a dipole (unnormalized) differential of the third kind [|3[ |29| with the 
residue resp ± dS = ^i§c{P±)dSo = ±iT at the first-order poles at P + and P_ on the curve 
1 6|). Let us introduce by now standard notation for the periods of the differential dSo 



a k = 6 dS , af = i dS . (4.14) 

Ja k J/34 
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One can verify from (|4.13|) , (|4.6| ) and (|4.4j ) that a k and a® take correspondingly real and pure 
imaginary values for the polynomial solutions. Using ( |4.9|) one obtains that the sum of all 
a— periods is given by the residue of dS at infinity 

ai + a 2 + ... + ajv-i = — 27U resp + dS = 2nT , (4-15) 

and only N — 2 periods are linear independent. The periods are functions of T and N — 2 moduli 
u k of the curve ( |4.6|) 

a k = a k (T;m, ...,u N -2) , k = 1,2, N - 2 (4.16) 

and a® has a similar dependence. 

Substituting the explicit expression ( |2.11| ) for the transfer matrix into ( |4.13| ) we can expand 



the meromorphic differential dSo over the basis of differentials of the first and the third kind on 
the curve Tn, oj k and dfl, respectively, as 



dS = T 

where 



N-2 

-2dVL + ( k + 2 ) u kduj k 

k=l 



dSl = T— , du k = T 1+k — , (k = l,...,N-2). (4.17) 

y y 

Here, additional powers of T were included to ensure invariance of the differentials under trans- 
formation ( [4.10| ), which acts as T — > XT and leaves u k invariant. The differentials ( |4.17| ) depend 



only on the moduli u k but not on T. Taking their linear combinations one can construct the 
canonical basis of the holomorphic 1-differentials du k with k = 1, N — 2 

N-2 . 

du k = Y u kj{u)dwj , f duj k = 2ixb jk (4.18) 

and the normalized differential of the third kind, dCl, with the residue resp ± dCl = ±1 at the 
first-order poles P± 

N-2 . 

dh = 2idfl + i Y Vj(u)dujj , f dQ = . (4.19) 

Here, U k j and Vj are real coefficients depending on the moduli of the curve Tn- Then, for given 
values of the periods ( |4.14| ) the differential dSo can be expanded over the canonical basis as 

N-2 

dS = iTdh + V — dCj k . (4.20) 
k=i 2n 

It is interesting to note that similar differential appears in the Seiberg-Witten description of the 
low-energy effective action of N = 2 SUSY Yang-Mills theory . The differential ( |4.2U| ) defines 
the asymptotic solution ( [3.30| ) of the Baxter equation, which being substituted into (|3.25|) , ( p.5| ) 
and C P . 1 3j ) determines the spectrum of the N Reggeon compound states. 

4.2. Hamiltonian flows 

The phase space for the system of N Reggeons is given by the direct product of the cycles a 3 - 
(j = 1,2, ...,N — 1) on the Riemann surface times the center-of-mass motion. The set of 
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points Qi, Qn-i on the Riemann surface situated one each on the a— cycles corresponds to 
the real values of the canonical Reggeon coordinates (xj,pj) and provides the coordinates on the 
level surface q k = const. Let us consider the Hamiltonian flow of the Reggeons on the curve T^, 
generated by the hamiltonians q k with the canonical Poisson bracket in the separated variables 
(xjiPj) as 



{Xj,q k } 



dq k 
d Pj 



with r k being the corresponding "times" . To calculate the r.h.s. we use the properties 
73]) of the separated variables to get 



2(coshp i - l)xf = q 2 x^- 2 + q 3 x^ :i + ... + q N 



N-2 



JV-3 



(4.21) 
and 

(4.22) 



with j = 1, N — 1 and invert these relations to find the expressions for the hamiltonians q k . 
It is convenient to introduce (N — 1) x (N — 1) matrix W inverse to the Vandermonde matrix 



x 



N-k 
3 



where 1 < j < N - 1, 2 < k < N, m(x 







-t k -i(x) 



(4.23) 



n 



N-l, 

i=i 



x 



^l<jl<:.<j k <N-l X h--- X 3k 



m'(xj) dxj 
- Xj) and t) 

symmetric polynomials. Then, one obtains from (|4.22| ) the expressions for the hamiltonians 

(4.24) 



q k 



N-l 

2 E W kj (x)x? (coshp, - 1) 

3=1 



which have a nonpolynomial kinetical part and which are very similar to analogous Hamiltonians 
for the classical Toda chain of iV interacting particles [p3| , |42j. The hamiltonians (|4.24j ) form 
a commutative family with respect to the Poisson bracket (|4.21| ). Substitution of ([4.24|) into 
( f4.21| ) yields the equations of motion for the Reggeon coordinates 



dxj 



2Wkj(x)xf sinhpj = W kj (x)JA 2 (xj) - Ax 



(4.25) 



where fl4.5|) and ( [4.3|) were used. These relations are well-known in the theory of finite-gap solu- 
tions of the KP/Toda systems as equations for zeros (or poles) of the Baker- Akhiezer function 

The integration of the evolution equations ( ^4.25| ) can be easily performed by the Abel map 
2^, |43|, |30fl . Indeed, multiplying the both sides of (|4.25| ) by x^~ l and taking into account that 



according to the definition ( f4.23| ), J2j W k jX 

dr k 



N-l 
3 



8 jh we get 



N-l 

E 

3=1 



A 2 (xj 



Axf 



(4.26) 



where k = 2, 
we obtain 



N. Comparing ( |4.26|) with the definition ( |4.17|) of the differentials duj k and Q 



t 2 = t 2 T 



N-l 

E 

3=1 



u 



r k T K 



k-l 



N-l 

E 

3=1 



du 



k -2 



(4.27) 



where k = 3, N and the points Qj lie on the Riemann surface one each over the a— cycles. 
Finally, we introduce a new set of the Reggeon coordinates 



(Qi, Qn-i) ~* (i?, <fi, fN-2) 



(4.28) 
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where 



r) 



N-l 



iV - 1 /"Qj 
XI / 



(4.29) 



and = 1, N — 2. Using the definitions (|4.27|) , (|4.18|) and ( |4.19| ) we find that the equations 
of motion ( 4.26 ) are trivially integrated in these variables 



N-2 

= 2t 2 + ]T Vi(u)t j+2 , 

.7 = 1 



7V-2 

V^fe = X U kj{u)t j+2 



(4.30) 



with {7 fcj - and Vj being t— independent functions of the moduli. 

The following remarks are in order. The number of points Qj in ( |4.28| ) is equal to the 
number of allowed bands, N — 1 = g + 1, and it does not match the genus of the Riemann 
surface T^. The vector with the coordinates <p = (ipi, ifN-z) depends on the integration path 
entering into (|4.29|) and it defines the point on the Jacobian <J(Tn) of the complex curve 



At the same time, the variable i? is related to the differential dCl and it starts to play a special 
role. Then, in new coordinates, i? and (<pi, fN-2) £ <J(Tn), the relations (|4.30|) describe 
a "fast" winding of the Reggeons around the Jacobi torus of the Riemann surface and a 



"slow" periodic motion in 1?. The corresponding periods can be evaluated using (|4.29|) , ( |4.18|) 
and C gTTSp as 



jk , 



dip k 



and 



d-d = 



dd 



ojv-i 
N-< 

*X 



iV-2 

X 



-2tt 



df2 — 2n resp, dfi = — 2n 



where j, k = 1, iV - 
are the hamiltonians ^3, 
in time t 2 . 



2. We recall that the generators of "fast" motion in times t 



3- 



t 



N 



q^-, while the hamiltonian q 2 generates "slow" drift of the system 



4.3. Action-angle variables 

The periods aj and the coordinates tpj have a simple interpretation in terms of the action-angle 
variables for the N Reggeon state. The differential dSo becomes the generating function of the 
canonical transformation from the separated coordinates (xj,pj) to the action-angle variables 
(<pj, Jj). Namely, the action variables are defined as 



Jk = — / dxp(x) = — i dS = 

2n Ja h 2lT Ja k 



Ok_ 

2tt 



' <*k 4" J C*k 

and the corresponding angles are given by 

dS r*k ^ cQ 



" r x k x rwk 

S = X / dxp(x)= X / dS 



a j ' 

UJk k=i " k=i 



where x^ belongs to the k— th allowed band. Substituting the canonical form ( 4.20|) of the 
differential dSo into this relation and taking into account that T = ^2^=1 Jk due to ( f4.15| ), one 
finds 

<fik = Vk - $ , 4>n-i = , f d(j)j> = 2n5 jr , (4.31) 
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where — 1, N — 1 and the last identity follows from ( |4.18|) and (|4.19| ). The angles 



>k 



describe the winding of the Reggeons around a— cycles on the Riemann surface and the 
corresponding basic oscillation frequencies are defined as 



e 



The frequencies 0^ do not depend on the times Tj and they can be easily evaluated from ( |4.31|) 



and ( f4.30| ) in terms of the coefficients U and V entering into definition of normalized differentials 
( gig ) and (14T9D Q. 



Summarizing our consideration of the dynamics of the N Reggeon system governed by the 
leading term in the asymptotic expansion of the Baxter equation solution (|3.30|) , we found that 
as a result of the composition of the maps 

fa,^) -> (xj,Pj) -> (Qj) e r N -+0,(<p x ,...,<p N -. 2 ) e J(T N ) 

the Reggeons have linear trajectories ( |4.30| ) on the Jacobian of the Riemann surface T^. Per- 
forming inverse transformations one can construct the explicit expressions for the holomorphic 
Reggeon coordinates by means of the Riemann theta functions corresponding to Tn E3L BUI. The 



resulting expressions are similar to the soliton solutions of the KP/Toda hierarchy [pO] , [29| , |43 | 
and in the center-of-mass frame of the N Reggeon compound state, P = 0, they can be repre- 
sented as 

z k = z k {t 2 , t 3 , t N ) = exp (2zt 2 ) • V k (t 3 U^ + ... + t N U^- 2 \ u h .., u N _ 2 ) . (4.32) 

Here, ^^ifij •••> <Pn-2] Wi, •••> w.w-2) is a 2vr periodic function of the variables y9j given by (|4.30| ). 
It depends on the parameters u k defined in ( |4.11| ) and on the (iV — 2)— dimensional wave vectors 
= (Uxk, ■ ■, Ujsr-2,k) built from the matrix U. In analogy with the KP/Toda solitons []30 



f4~3| , |29| , the wave vectors can be expressed in terms of the (3— periods of the differentials of the 
second kind, dti k , on the curve T^, normalized as E3 



dflk = , dCl k 



where £ = T/x — > is a local parameter on the curve at the vicinity of infinity P + on the 
upper sheet. Using the well-known property [f||, that the (3— periods of the differentials dCt k 



are related to the behaviour of the holomorphic differential dCj k near infinity P + as 

j dn k = i^-'/iCO u > ~° 



and substituting ([4.18 ) and ( |4.17|) into this relation, one can obtain the following expressions 



for the coefficients U k j = lljf^ 



Ui 1} =2i£ dtii , Ui 2) + ^Ul 1} =Ai£ dCl 2 , ... . (4.33) 

Jp k 2 Jp k 



7 The same expression can be found using the fact that the transition from separated coordinates (xj,pj) 
to the action-angle variables (J k ,y>k) is the canonical transformation and the evolution of if is described by the 
Hamiltonian equation 

— = {v k , qj } = — = 2«— = e kj , 



where periods were defined in ( 4.14 ) and ( 4.16 ) 
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To understand the origin of the first factor in ( |4.32|) involving the time t 2 = t 2 T, we observe 



that (|4.32|) gives an exact solution to the hierarchy of the conservation laws 

Let us consider the evolution of z k in time r 2 . Replacing in (|2.14j) ttj = —irjd Zj by their classical 
analogs, one finds that the Hamiltonian q 2 is related to the holomorphic part of the total angular 
momentum of the N Reggeon system in the center-of-mass rest frame, Y^j Kj — P — 0, as 

<?2 = -T 2 = z k^kj ■ 

Therefore, the evolution of the Reggeon holomorphic coordinates in time t 2 becomes trivial 

d T2 z k = {z k ,q 2 } = 2iTz k =>> z k ~ exp (2ir 2 T) 
and it corresponds to the rotation of the N Reggeon system on the 2-dimensional plane of the 



impact parameters ( |2.2|) with the angular momentum 2T. Then, using the definitions (|4.12j ), 



( |2.14 ) and ( |2.15| ) one obtains that in the large h limit the total angular momentum of the soliton 



wave is equal by the Lorentz spin m of the N Reggeon state 

1 



2T = 2y/h(h - 1) =2[h-^) + O^ 1 ) =m. 



Thus, the first factor in (|4.32 ) describes the rotation of the N Reggeon system around its center- 



of-mass in the 2-dimensional impact parameter space. For N = 2 Reggeon state, the BFKL 
Pomeron, this becomes the only mode of the Reggeon motion. For N > 3 the quasiclassical 
dynamics of the N Reggeon compound states becomes much more interesting due to appearance 
of the soliton excitations ( |4.32j ). The space of the parameters u±, un-2 of the N Reggeon 
soliton waves coincides with the moduli space M. of the hyperelliptic curve T^- The coordinates 
on M. are determined by the values of the quantum numbers q 2 , qN- 



4.4. Singularities on the moduli space 



BTL EH 



indicate that for fixed q 2 the possible 

2)— dimensional 



The numerical solutions of the Baxter equation 

values of the quantum numbers qs, q^ occupy a compact region in (N 
space with the boundary defined by the hypersurface (|3.22|) . Let us show that the hypersurface 
( |3.22| ) can be identified from the analysis of singularities on the moduli space M. of the N 

] that the curve Tn becomes singular only when 



Reggeon soliton waves. It is well-known |43 
two branch points merge, Oj 



Oj+i (see Fig. 1). Since the positions of tjj depend on the values 
of the quantum numbers, the latter relation implies certain conditions on q k . 
For N = 2 Reggeon states, one finds the branch points of the curve ( |4.7| ) as 



-01 



T 

a 2 = — . 
2 



Since T 2 = — q 2 > for the polynomial solutions, the values of 01 and a 2 are always real and 
different^. Thus, the curve T 2 does not have any singularities and this is in agreement with the 

8 We recall that the band [a%, 02] defines the interval of classical motion in the separated coordinate X\. 
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fact || |26| that the exact solution of the Baxter equation for iV = 2 is well defined and the 
asymptotic expansion of the energy (|3.21|) is Borel summable for N = 2. 

For N = 3 Reggeon states, or the QCD Odderon, the curve ( |4.8| ) has four branching points^ 



4r? + g 2 o- i + g 3 = 0, (j = 1,2,3) 



04 = -- 



03 

fa 



To classify all possible solutions we introduce the "effective" discriminant A of the curve r 3 as 
A = 16a 4 Vi - <r 2 ) V2 - ^3)Vs - oi? = q\ u\{l - 21u\) , 



where u\ was defined in ( [4. 1 1| ) . Then, for A > all roots are real; for A < two roots are 
real and two remaining roots are complex conjugated to each other and, finally, for A = two 
roots coincide and the curve r 3 becomes singular. 

Therefore, in order to be able to construct two real intervals corresponding to the classical 
motion of the Reggeons in the separated coordinates x\ and X2 one has to require A > 0. If one 
assumes that the moduli space M. is the complex u\— plane, then the singularities are located 
at three points 

1 , and -^L (4.34) 



ut nt 



'27 V27 

and the polynomial solutions of the Baxter equation for N = 3 Reggeon states correspond to 
the real values of U\ shown in Fig. 2 such that 



< u\ < 



1 

27 



(4.35) 



For A = — 1 the property ( |4.10| ) of the curve T 3 leads to the symmetry on the moduli space M. 
under u\ <-> —u\. 




Figure 2: The possible values of U\ for iV = 3 Reggeon states corresponding to the polynomial 
solutions of the Baxter equation. The dots indicate the positions of singularities. 



The relations ( |4.34| ) and (^4.35|) are in complete agreement with the numerical results [^1], |26 



and with asymptotic expansions of the solutions of the Baxter equation, (|3.23|) . A similar 



analysis can be carried out for higher iV > 3 Reggeon states p6 |. 



3 Here we do not require the ordering o\ < 01 < 03 < 0-4. 
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Any point on the moduli space M. corresponds to a certain Riemann surface and it 
describes the compound state of N Reggeons with the quantum numbers defined by (|4.12j ). 
The energy e n of these states becomes the function on the moduli space and it is naturally 
to expect that singularities on the moduli space M. control the analytical properties of the 
functions for iV > 3. They are responsible for the appearance of Borel singularities in the 
asymptotic expansion of the energy (|3.21|) for iV > 3. 



5. Quantization conditions 

Let us discuss the quantization conditions for eigenvalues of the integrals of motion q 2 , 



q jy. The quantization of q 2 follows from the definition ( [2.141 ) after one takes into account that 
possible values of the conformal weight are given by (|2.15| ). The polynomial solutions of the 
Baxter equation correspond to the special values ( |3.14| ) of quantized q 2 = —h(h — l). In this case 
one has to establish the quantization conditions for the remaining charges qs, qN and then 



try to generalize them from integer conformal weight h to all possible complex values ( [2.15 ). 



Inverting the dependence (|4.16| ) and using (|4.15| ) one can obtain that the quantized values 



of the moduli ( |4.11| ), or equivalently qs, qN, are determined by q 2 and by iV — 2 linear 
independent periods a k of the differential dSo defined in ([4.14]) 

u k = u k (q 2 ;ai,a 2 , ...,a N - 2 ) , k = 1, N - 2 . (5.1) 

Although q 2 and a k depend separately on the parameter r], this dependence is cancelled inside 
the function u k due to invariance of the moduli under transformation Q4.10 ). 



The periods a k have a simple interpretation in terms of the roots Xj of the solutions Q3.15| ) 



of the Baxter equation. We recall that the solutions Q(x) define the wave function of the N 
Reggeon state in the separated coordinates ( |3.5| ). Therefore, the roots Aj of Q(x) being the 
zeros of the Reggeon wave function should belong to the intervals of the classical motion of 
Reggeons, that is to the N — 1 allowed bands on the Riemann surface T^- Let us denote by n k 
the number of roots of Q(x) (including the N— time degenerate root at x — 0) which belong to 
the k— th allowed band, Xj G [o~ 2k -i,<J 2k ], 



N-l 



n k >0, J2 n k = h ( 5 - 2 ) 



k=l 



and rij > N if the root x = belongs the j— interval <J 2 j-i < < a 2 j. Then, it follows from 
( |3.28| ) that the meromorphic differential dS has first-order poles on the Riemann surface r^v at 
x = Xj and its periods around a— cycles count the number of roots 

f dS = 2nr]n k . 
Substituting the expansion ( |3.27| ) of dS in powers of 77 one gets 

£ dS = 2tt V (n k + -) + 0{rj 2 ) , (5.3) 



where 0(rf) term takes into account the contribution of higher terms in the expansion ( 3.27|) . 



The relations (|5.3| ) take the form of Bohr-Sommerfeld quantization conditions for the Reggeon 
wave function ( |3.30| ). 
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Substituting flOl) into ( |5.1|) we find that the moduli of the Riemann surface, or equivalently 
the integrals of motion g 3 , q^, become quantized and, in accordance with our expectations 
( |3.17 ), their values are parameterized by q2 and by the set of positive integer numbers n k . 



However, trying to find the dependence of moduli on n k from (|5.1|j and ( |5.3|) , one has to take 
into account that the periods a k themselves are functions of the moduli due to 0(rj 2 ) term in 
([5.31), which does not depend on integers n k and has the general form rf(— q 2 ) 1 ^ 2 H{{u k }). Here, 
H is a complicated function of the moduli and q 2 is needed to restore the scaling dimension 



of a k in rj under transformation (|4.10| ). Then, using independence of u k on rj one can put 
rj = (— q-if 1 ^ 2 in (|5.1| ) and represent the general solution of ( |5.1|) and (|5.3j ) as 



TbA 1 1 



U k = U k (d 1 ,8 2 , ■■■^N-2\q 2 l ) , 8j = - 3 ,i/2 • ( 5 - 4 ) 

(-92) 

One possibility to define the function u k is to find its asymptotic expansion in inverse powers 
of the conformal weight using the series ( |3.21|) , 

u k = Yl c k( n ) h ~ k ( 5 - 5 ) 

fc>0 

with q 2 = —h(h — 1) and rij = fixed. However, due to the presence of singularities on the 
moduli space A4, the series ( |5.5|) turns out to be non Borel summable. Let us now change the 
parameters of the expansion and expand the function ( |5.4j) in powers of q 2 l keeping 5j = fixed. 



As example, one uses the result [26] for the large h expansion of u\ to 0(h ) order for N = 3 



Reggeon states to convert it into the following form 



u 1 = Y,< ) (S)q 2 k , (5.6) 

fc=0 



where the leading term is given by 
^(l-35 + 25 2 - 

Tt V 9 81 243 2187 6561 



4 0) (6) = 4= ( 1 " 3 S + 2 6 2 - -6 3 + -5 4 + + + + 0(5 8 ) ) . (5.7) 



We notice that the expansion ( |5.6| ) goes over integer powers of l/q 2 = 0(h~ 2 ), while the original 
series (|5.5| ) had a much bigger parameter of the expansion, 1/h. 

Let us consider now the asymptotic approximation to solution of the Baxter equation given 
by the expression ( |3.30|) , in which we neglect all nonleading 0(rj) corrections to the exponent. 
In this limit, there are no 0(rf) corrections to the periods a k in ( |5.3| ) and the expressions ( p.l[ ) 
and ( |5.4|) for quantized moduli look like 

u[ 0) =u k (q 2 ;ai,a2,...,a N -2)\ ak=2wv f nk+ i\ = u k (S h 5 2 , 6 N - 2 ; 0) , (5.8) 



where the parameters S k were defined in (]5. 

To obtain all possible values of quantum numbers q^, q^ one has to evaluate the moduli 
for different sets of integers n\, n N _ 2 and substitute them into (|4.12j) . Let us consider 



as a definition of a continuous function of real positive q 2 , which for q 2 = —h(h — 1) and 



h positive integer gives quantized q k . Then, for different possible sets ( |5.2j ) of integers n k the 
functions ( |5.8| ) define the family of curves on the moduli space A4. Each curve describes the 
flow on the moduli space Ai in the "slow" time q 2 and it has a distinguished property that the 
values of the periods a k are preserved. As example [p6~| , the flow of quantized ^3 as a function 
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zoo- 




Figure 3: One of the curves describing the flow of quantized q% for N = 3 Reggeon states. 
The circles represent the results of numerical calculations. Three dashed lines correspond to the 
singularities on the moduli space (|4.34|) . 



of q2 for N = 3 Reggeon states and fixed value of integer n\ = 4 is shown by solid line in Fig. 3. 
This curve induces the flow on the moduli space indicated by the solid line in Fig. 2. 

For N = 3 Reggeon states, is equal to the leading term of the expansion (|5.6|) given by 



Comparing (|5.7|) with ( 4.34J) we find that the expression (|5.7|) provides a weak coupling 
expansion of the moduli in 5 around one of the singularities on the moduli space in Fig. 2. To 
approach two remaining singular points one has to develop the strong coupling expansion of the 
moduli. 

The Seiberg-Witten formalism gives us a powerful method of calculating the moduli (|5.8| ), 
based on the Whitham equations and on the remarkable property of duality between strong and 



weak coupling expansions of the moduli Uk in parameters 8j [44]. In application to the N = 3 
Reggeon states, the duality originates from the property @ that the quantity (ai(ui), af(ui)), 
built from the periods ( |4.14|) of the differential dSo on the curve T3, has well-known monodromies 
around the singular points (|4.34|) . The monodromies are given by 2 x 2 matrices, which belong 
to the subgroup T(2) of the SL(2, Z) group consisting of the matrices congruent to 1 modulo 
2.0 The symmetry of the spectral curve T 3 under n x — > —U\ leads to the following property of 
the moduli JI3| 

«i(l - 5) = -mi(<5) , 



which together with (|5.7|) allows us to identify the values of S corresponding to the singularities 
( [4.34|) on the moduli space as 

5 crit = 1 , - and , 
' 2 

respectively. Having expressions for the monodromy of (ai(wi), a^iui)) one can determine the 

10 This explains the observation made in |27j] that the Schrodinger equation for the N — 3 Reggeon state, QCD 
Odderon, obeys a new modular symmetry with respect to T(2). 
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asymptotic behaviour of the moduli around these points in the following form 



"44 



S _ i „ -V3^ x + -L^ , <f-^~^ui(lnui-l), 5 ~ -V3^ - . (5.9) 



These relations describe the flow of the quantum numbers shown in Fig. 3 and they can be also 
derived from the Whitham equations. 



5.1. Whitham equations 



Let us show that the flow of the moduli (|5.8j ) in the "slow" time - variable q 2 is governed by 
the Whitham equations JTT|, |9|. We recall that the function ( |5.1| ) is inverse to ( |4.16| ) and the 
dependence of on q 2 in ( |5.8|) can be found from the condition that the periods should be 
invariant under variations of q 2 . Using the definition of the periods ( 4.14J ), this condition can 
be expressed as follows 



8a 



k = ST f d T dS + Suj f d Uj dS = . 



(5.10) 



To calculate the external derivatives one considers the variation of the differential dSo with 
respect to Reggeon quantum numbers. As follows from the definitions Q4.13Q and 



„„ , Scu dx 
oabo = ax — = dA(x) — 
u y 



] 2 X 



N-2 



+ 5q 3 x N ~ 3 + ... + 5q 



N 



dx 



where we put rj — 1 for simplicity. Finally, one gets 

ddS „ , ddS 



d Uj 



Tdoo 



3 i 



dT 



T- l dS . 



where holomorphic differentials dujk were defined in ( |4.17|) . Here, the first relation states that 
the variation of the differential dSo with respect to moduli is proportional to the holomorphic 
differential defined on the curve T^- Being applied to ( 5.10|) , this remarkable property of the 
differential dSo leads to the Whitham equations for the moduli 



JV-2 



dui 



"E^TH dcoj+f dS = 

i=1 Ja k Ja k 



(5.11) 



Calculating the a— periods of the differential duk from ( |4.18D as § douj = 2n(U 1 ) J -fc and taking 
into account ( |5.3| ), we obtain the following system of equations 



T 2 du i 

dr 



N-2 



k=\ 



1 



J2 \ n k + 2j u kj(u) , 



(5.12) 



where j = 1, N — 2 and Ukj are functions of the moduli u\, un-2- The matrix elements 
of Ukj define the wave vectors of the soliton waves Q4.32 ) and one can rewrite (|5.12 ) as 

where notation was introduced for the vector v k = + | and (v, U^) = v^Ukj- 
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For N = 3 Reggeon states, the system (|5.12|) is reduced to an ordinary differential equation 
for the moduli u\ = Ui(8) with 6 = (n+ |)/T 



(5.13) 



where U = U {1) was defined in ( [4.33|) . It allows us to determine exactly the function entering 
into Q5.7 ) and study its properties at the vicinity of singularities ( |3.22| ) on the moduli space Ai. 
One can show |Q that the Whitham equation ( |5.13| ) is in agreement with the weak coupling 
expansion (|5.7|) and with the asymptotics Q5.9|) . 



5.2. Boundary conditions 



To solve the differential equations ( |5.12|) and ( |5.13|) one has to supplement them by an appro- 



priate boundary conditions on u^. These conditions are provided by the asymptotic behaviour 
of Uk as q 2 — > — oo, which follows from the large h expansion of quantum numbers (|3.21| ). 



Let us start with N = 3 Reggeon states and examine the flow of quantized q 3 shown by 
solid line in Fig. 3. For — q 2 = h(h — 1) — > oo the quantized q 3 behave as g 3 ~ (— g 2 ) 3 ^ 2 /v / 27 
leading to the asymptotic expression for the moduli u\ = -7=, which we identify as one of the 
singularities on the moduli space M.. Thus, the singularity on the moduli space, 

Ui(<J) = -^=, for 5 = (5.14) 



becomes the starting point for the Whitham evolution ( |5.13| ). As — q 2 decreases toward the 
origin, u\ passes the second singularity (|3.22| ) at U\ = 0, and for smallest — q 2 , corresponding 
to the boundary of the polynomial solutions, u\ approaches the third singularity on the moduli 
space at U\ = — 77^- We conclude that for polynomial solutions of the N = 3 Baxter equation, 



the Whitham equation ( 5.13 ) describe the trajectory on the moduli space (see Fig. 2), which 
goes along the real axis from ^7= to — The Whitham evolution of the moduli corresponding 
to higher A^ > 3 Reggeon states follows a similar pattern. For — q 2 — > 00 the quantum numbers 
qk approach the values belonging to the critical hypersurface ( |3.22|) . The same values coincide 



with the positions of singularities on the moduli space M., Uk ~ u^ ng as — q 2 — > 00. Since 
the possible values of the moduli corresponding to the polynomial solutions occupy a compact 
region on the moduli space with the boundary being the singularities, the flow of Uk will start 
at one singular point as —q 2 — > 00 and will finish at another singular point on M. for — q 2 taking 
the smallest value allowed for the polynomial solutions of the Baxter equation. 

As was shown in Sect. 4.2, the "time" variable in the Whitham evolution, T, is related to 
the total angular momentum associated with the rotation of the A^ Reggeon system around 
its center-of-mass in the impact parameter space. The Whitham equations ( 5.12|) and (|5.13 ) 



describe the adiabatic perturbation of the moduli of the curve which enter as parameters 
into the soliton solutions (|4.32|) for the A^ Reggeon states. Thus, in the leading nonlinear WKB 
approximation, (|3.30|) , the A^ Reggeon states corresponding to the integer positive conformal 
weight h can be considered as modulated soliton waves. 

We recall that till now we worked on the subspace of the polynomial solutions of the Baxter 
equation. The Whitham equations ( |5.12j ) and ( p.l^| ) offer a natural way of analytical continua- 
tion of the obtained expressions beyond this subclass. Let us consider for simplicity the A^ = 3 
Reggeon states. The evolution of q 3 is described by a smooth function of q 2 = —h(h — 1) (see 
Fig. 3), which for h in Q3.14Q gives the values of quantized q 3 . We observe that, first, the same 
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function allows us to define formally q% corresponding to half-integer h in ( [2.15| ) and, second, the 
evolution of ^3 inevitably leads to the region of small and positive q?,, where polynomial solutions 
do not exist. Applying the Whitham equations ( |5.12| ) and ( |5.13| ) in these two cases one assumes 
that the above interpretation of the N Reggeon states as modulated soliton waves holds for any 
integer and half-integer conformal weight h. One may also try to apply the Whitham equations 
for an arbitrary complex quantized conformal weights ( |2.15|) , but the important difference with 
real h is that one has to find now the initial conditions similar to ( |5.14| ). 



6. Conclusions 

The Regge asymptotics of hadronic scattering amplitudes in high-energy QCD are controlled by 
the color-singlet compound states of Reggeons. Reggeons appear as a new collective degrees of 
freedom of QCD in the Regge limit and their dynamics in four dimensions is described by the 
effective (1 + 1) -dimensional Hamiltonian, which exhibits remarkable properties of integrability. 
The system of N interacting Reggeons in the multi-color limit resembles very much the 1- 
dimensional Heisenberg spin s = chain with iV sites. It has enough number of conserved 
charges qk to be completely integrable. To diagonalize the N Reggeon hamiltonian and calculate 
the spectrum of the Reggeon compound states, QCD Pomerons and Odderons, we defined a 
new set of Reggeon coordinates, in which N coupled Schrodinger equations for eigenvalues of 
the conserved charges qk become separated and are replaced by the Baxter equation (|3.12j ). 
The solutions of the Baxter equation, depending on the set of quantum numbers qk, define 
the energy and the wave function of the N Reggeon state in the separated variables. It is the 
Baxter equation that summarizes the QCD dynamics of Reggeons in the separated coordinates 
and whose nonlinear WKB expansion gives rise to the integrable structures well-known from 
the finite-gap solutions of the soliton equations and their Whitham deformations. 

The leading term of the nonlinear WKB expansion of the polynomial solutions of the Baxter 
equation defines the hyperelliptic Riemann surface IV as the level surface of the integrals 
of motion qk = const, and the meromorphic 1-differential dSo on it. The conserved charges 
qk generate the Hamiltonian flows of Reggeons on and the exact solution of the arising 
hierarchy of the evolution equations is given by the Reggeon soliton wave. The moduli Uk of 
for N > 3 depend on the quantum numbers of the Reggeon states and enter as parameters into 
the soliton solutions. The possible values of q%, qx are quantized and they determine the 
family of curves on the moduli space Ai. Each curve describes the flow on Ai, which is governed 
by the Whitham equations. These equations describe the adiabatic perturbation of the Reggeon 
soliton waves and the properties of their solutions will be considered in the forthcoming paper 

II- 

We would like to mention in conclusion that had we perform a similar nonlinear WKB 
analysis of the Baxter equation for Toda chain ( |3.24| ) in the special case % = 1, we could 
reproduce the main ingredients of the Seiberg-Witten solution || of the effective action of 
N = 2 supersymmetric QCD. It remains unclear, however, what are a new collective variables 
in N = 2 SUSY QCD, which play a role similar to the separated coordinates of Reggeons in 
high-energy QCD. 
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